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Abstract 

We study exact multi-soliton solutions of integrable hierarchies on noncommutative 
space-times which are represented in terms of quasi-determinants of Wronski matrices by 
Etingof, Gelfand and Retakh. We analyze the asymptotic behavior of the multi-soliton 
solutions and found that the asymptotic configurations in soliton scattering process can 
be all the same as commutative ones, that is, the configuration of iV-soliton solution has 
N isolated localized energy densities and the each solitary wave-packet preserves its shape 
and velocity in the scattering process. The phase shifts are also the same as commutative 
ones. Furthermore noncommutative toroidal Gelfand-Dickey hierarchy is introduced and 
the exact multi-soliton solutions are given. 
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1 Introduction 



Extension of integrable systems and soliton theories to non-commutative (NC) space- 
times § have been studied by many authors for the last couple of years and various kind 
of integrable-like properties have been revealed [U [2]. This is partially motivated by 
recent developments of NC gauge theories on D-branes. In the NC gauge theories, NC 
extension corresponds to introduction of background magnetic fields and NC solitons are, 
in some situations, just lower-dimensional D-branes themselves. Hence exact analysis of 
NC solitons just leads to that of D-branes and various applications to D-brane dynamics 
have been successful [3j. In this sense, NC solitons plays important roles in NC gauge 
theories. 

Most of NC integrable equations such as NC KdV equations apparently belong not 
to gauge theories but to scalar theories. However now, it is proved that they can be 
derived from NC anti-self-dual (ASD) Yang- Mills equations by reduction [4J, which is 
first conjectured explicitly by the author and K. Toda [5]. (Original commutative one 
is proposed by R. Ward [6] and hence this conjecture is sometimes called NC Ward's 
conjecture.) Therefore analysis of exact soliton solutions of NC integrable equations could 
be applied to the corresponding physical situations in the framework of N=2 string theory 

mum- 

Furthermore, some soliton equations describe real phenomena such as shallow water 
waves in fluid dynamics, optics and so on. If noncommutativity in space-time affects 
soliton dynamics, then we can check whether our universe is noncommutative or not by 
comparing experimental results and estimate the strength or the upper bound of the 
noncommut at i vity 

Hence, construction and analysis of exact multi-soliton solutions are worth studying 
from various viewpoints of integrable systems, string theory, and perhaps detection of 
noncommutativity in our universe. 

Exact multi-soliton solutions of noncommutative KP hierarchy are constructed by 
Etingof, Gelfand and Retakh in 1997 [10J, where quasi-determinants play crucial roles. 
(For other applications of quasi-determinants to noncommutative integrable systems, see 
e.g. [m d2l EH El]-) However, their discussion is general and explicit analysis of the 
behavior of their soliton solutions has not yet been done. Paniak also constructs multi- 
soliton solutions of NC KP and KdV equations (not hierarchies) and studies the scattering 

2 In the present paper, the word "NC" always refers to generalization to noncommutative spaces, not 
to non-abelian and so on. 
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process [15] a However, the discussion about the soliton dynamics is mainly focused on 
two-soliton scatterings. 

In this paper, we study exact multi-soliton solutions of NC integrable hierarchies in 
terms of quasi-determinants of Wronski matrices, which is developed by Etingof, Gelfand 
and Retakh. We analyze the asymptotic behavior of the multi-soliton solutions and found 
that the asymptotic configurations can be real-valued though NC fields take complex 
values in general. The behavior in soliton scatterings is all the same as commutative 
ones, that is, the iV-soliton solution has N isolated localized energy densities and the 
each wave-packet preserve its shape and velocity in the scattering process. The phase 
shift is also the same as commutative one. 

This paper is organized as follows. In section 2, we make a brief introduction to NC 
field theory in star-product formalism. In section 3 and 4, we review definition and some 
properties of quasi-determinants and their applications to construction of multi-soliton 
solutions of NC integrable hierarchy in star-product formalism. In the end of section 4, 
we introduce NC toroidal Gelfand-Dickey (GD) hierarchy and give exact multi-soliton 
solutions which are new. In section 5, we discuss the asymptotic behavior of them in 
detail. Section 6 is devoted to conclusion and discussion. 

2 NC Field Theory in the star-product formalism 

NC spaces are defined by the noncommutativity of the coordinates: 

[x i ,x>]=i6 ii , (2.1) 

where the constant d v is called the NC parameter. If the coordinates are real, NC pa- 
rameters should be real. Because the rank of the NC parameter is even, dimension of NC 
space-times must be more than two. Hence in this paper, we deal not with integrable 
systems in (0 + l)-dimension such as the Painleve equation, but with ones in (1 + 1) 
or (2 + l)-dimension such as the KdV and KP equations. In (1 + l)-dimension, we can 
take only space-time noncommutativity as [t,x] = id. In (2 + l)-dimension, there are 
essentially two kind of choices of noncommutativity, that is, space-space noncommuta- 
tivity: [x, y] = id and space-time noncommutativity: [t, x] = id or [t, y] = id, where the 
coordinates (x, y) and t correspond to space and time coordinates, respectively. 

NC field theories are given by the replacement of ordinary products in the commu- 
tative field theories with the star-products and realized as deformed theories from the 

3 Dimakis and Miiller-Hoissen present perturbative corrections with respect to a noncommutative 
parameter in 2-soliton scatterings of the NC KdV equation [16] before the Paniak's work. 
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commutative ones. The star-product is defined for ordinary fields on flat spaces, explic- 
itly by 



[TT] . The ordering of fields in nonlinear terms are determined so that some structures such 
as gauge symmetries and Lax representations should be preserved. 

The star-product has associativity: f*(g*h) = (f*g)-kh, and reduces to the ordinary 
product in the commutative limit: — > 0. The modification of the product makes the 
ordinary spatial coordinate "noncommutative," that is, [a;*, a;- 7 ]* := x % -kx^ — x^-kx 1 = i9 l i . 

We note that the fields themselves take c-number values as usual and the differenti- 
ation and the integration for them are well-defined as usual. A nontrivial point is that 
NC field equations contain infinite number of derivatives in general. Hence the integra- 
bility of the equations are not so trivial as commutative cases, especially for space-time 
noncommut at i vity. 

3 Brief Review of Quasi-determinants 

In this section, we make a brief introduction of quasi-determinants introduced by Gelfand 
and Retakh [18, 19] and present a few properties of them which play important roles in 
the following sections. The detailed discussion is seen in e.g. [201 !2"Tj . Relation between 
quasi-determinants and NC symmetric functions is seen in e.g. |22j. 

Quasi-determinants are not just a generalization of usual commutative determinants 
but rather related to inverse matrices. From now on, we suppose existence of all the 
inverses. 

Let A = (ciij) be a. N x N matrix and B = (fry) be the inverse matrix of A, that is, 
A-kB = B-kA = 1. Here all products of matrix elements are supposed to be star-products, 
though the present discussion hold for more general situation where the matrix elements 
belong to a noncommutative ring. 

Quasi-determinants of A are defined formally as the inverse of the elements of B = A^ 1 : 






where d\ 



dj dx H and so on. This explicit representation is known as the Moyal product 




In the commutative limit, this is reduced to 




(3.2) 
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where A % i is the matrix obtained from A deleting the i-th row and the j'-th column. 

We can write down more explicit form of quasi-determinants. In order to see it, let us 
recall the following formula for a block-decomposed square matrix: 

(A B\ X _( {A - B * D- 1 * C)- 1 -A~ l *B*(D -C + A'^B)- 1 \ 

[C D ) ~ y —(D — C*A~ 1 *B)~ 1 *C*A~ 1 (D - C * A' 1 * B)- 1 J' 

where A and D are square matrices. We note that any matrix can be decomposed as 
a 2 x 2 matrix by block decomposition where one of the diagonal parts is 1 x 1. Then 
the above formula can be applied to the decomposed 2x2 matrix and an element of the 
inverse matrix is obtained. Hence quasi-determinants can be also given iteratively by: 



\A\ij = a i:j - a a' * ((A tJ ) ^i'j' + a 

= a ij - aii ' * d^lfi')' 1 * a 



3 3 ' 



(3.3) 



It is sometimes convenient to represent the quasi-determinant as follows: 

an • • • a>ij ■ • • Qin 



I n 



an 



"13 



a n i ■ ■ ■ a 



"j 



(3.4) 



Examples of quasi-determinants are, for a 1 x 1 matrix A = a 

\A\ = a, 



and for a 2 x 2 matrix A = (aij) 



\A\ 
\A 



ii 



21 



an 


Ol2 


On 


«22 


an 


Ol2 


021 


«22 



an — ai2 * a 2 2^ * a2i, \A 



— a2i — a22 * o 12 * an, 



12 



22 



an 


ai2 


021 


<*22 


an 




«21 


022 



a 12 — Oil * 021 1 * °22 ? 



— ^22 — 021 * a n * ai2, 



and for a 3 x 3 matrix A = (a^) 

i Qn | a 12 &i3 
W|ll = 021 022 023 
Q31 ^32 ^33 



an - (a 12 ,a 13 ) * 



a ll ~~ O12 * 



a i3 * 



&22 «23 



«32 033 

022 ^23 
«32 &33 



-1 



* a 2 i — ai2 * 



^22 023 
a 32 &33 



^22 0-23 



021 

a 3 i 



a 32 0-33 



-1 



^ a 2 i - a 13 * 



^22 ^23 
0-32 0-33 



*a 3 i 



★ a 3 i, 
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and so on. 



4 Exact Soliton Solutions of NC Integrable Hierar- 
chies 

In this section, we give exact multi-soliton solutions of several NC integrable hierarchies in 
terms of quasi-determinants. In the commutative case, determinants of Wronski matrices 
play crucial roles. In the NC case, these determinants are just replaced with the quasi- 
determinants. We review foundation of the NC KP hierarchy and the /-reduced hierarchies 
(so called NC GD hierarchies or NC /KdV hierarchies), and present the exact multi-soliton 
solutions of them developed by Etingof, Gelfand and Retakh |10j . Finally we extend their 
discussion to the NC toroidal GD hierarchy. 

An iV-th order pseudo-differential operator A is represented as follows 

A = a N d x + CLN-xd^' 1 H h a + a-\d~ l + a_ 2 <9~ 2 H , (4.1) 

where a, is a function of x associated with noncommutative associative products (here, 
the Moyal products). When the coefficient of the highest order a at equals to 1, we call it 
monic. Here we introduce useful symbols: 

A> r := d™ + a^idf- 1 + • • • + a r d r x , (4.2) 
A< r := A — A> r+1 = a r d r x + a r _i<9; _1 + ■ • • . (4.3) 

The action of a differential operator <9™ on a multiplicity operator / is formally defined 
as the following generalized Leibniz rule: 



{d x f)d n ~\ (4.4) 



where the binomial coefficient is given by 

n(n — 1) • • • (n — i + 1) 



n 



i ' 1) •• -1 



(4.5) 



We note that the definition of the binomial coefficient (14.51) is applicable to the case for 
negative n, which just define the action of negative power of differential operators. 

The composition of pseudo-differential operators is also well-defined and the total set of 
pseudo-differential operators forms an operator algebra. For a monic pseudo-differential 
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operator A, there exist the unique inverse A -1 and the unique m-th root A l l m which 
commute with A. (These proofs are all the same as commutative ones.) For more on 
pseudo-differential operators and Sato's theory, see e.g. j23j [2U [25l [26] . 

In order to define the NC KP hierarchy, let us introduce a Lax operator: 

L = d x + uid- 1 + u 3 d~ 2 + u 4 d~ 3 + ■■■, (4.6) 

where the coefficients Uk (k — 2, 3, . . .) are functions of infinite coordinates x := (xi, x 2 , . . .) 
with x\ = x: 

u k = u k (x 1 ,x 2 , • . •)• (4.7) 

The noncommutativity is introduced into the coordinates (xi,x 2 , . . .) as Eq. (12. f I) here. 
The NC KP hierarchy is defined in Sato's framework as 

d m L=[B m ,Ll, m = l,2,..., (4.8) 

where the action of d m := d/dx m on the pseudo-differential operator L should be in- 
terpreted to be coefficient- wise, that is, d m L := [<9 m ,L]* or d m d^. = 0. The differential 
operator B m is given by 

B m := ( L*---*L )> =: (L m )> . (4.9) 

m times 

The KP hierarchy gives rise to a set of infinite differential equations with respect to 
infinite kind of fields from the coefficients in Eq. (14. 8[) for a fixed m. Hence it contains 
huge amount of differential (evolution) equations for all m. The LHS of Eq. (14. 8 j) becomes 
d m Uk which shows a kind of flow in the x m direction. In the a^-flow equations, we can 
see that infinite kind of fields 1*3, W4, W5, ... are represented in terms of one kind of field 
u 2 [271 [28]. 

If we put the constraint (L')<_i = or equivalently L l = Bi on the NC KP hierarchy 
(14. 8p . we get a reduced NC KP hierarchy which is called the l-reduction of the NC KP 
hierarchy, or the NC I KdV hierarchy, or the Z-th NC G elf and- Dickey hierarchy. Especially, 
the 2-reduction of the NC KP hierarchy is just the NC KdV hierarchy. Explicit examples 
are seen in e.g. [28]. (See also [291 EH ED]-) 

Now we construct multi-soliton solutions of the NC KP hierarchy. Let us introduce 
the following functions, 

f s (x) = ei^' as) + a s ei& A \ (4.10) 
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where 



a) = x\a. + x 2 a 2 + x^a 3 + • • • , 
and a s , f3 s and a s are constants. Star exponential functions are defined by 

OO 1 



n=l 



n times 



An iV-soliton solution of the NC KP hierarchy ( 14.81) is given by [10J, 

L = <5> N *d x <5> N \ 

where 



./ 





■ ■ , In, f) \n+i,n+i, 


h 
fl 


h ■■■ 
fk ■■■ 


fs 
In 


AN-l) 

k (2v) 


AN-l) 

J2 
AN) 
J2 


AN-l) 
J N 

AN) 

Jn 


■ j fm) 


is given by 




( 


ft h 





f(N-l) 



fi N) 



fl f-2 



fm \ 
f 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



(4.15) 



. f(m-l) I 
J m > 

where /i, f 2l fm are functions of x and /' := df/dx, f":=d 2 f/dx 2 , f^ :=d m f/dx 
and so on. 

In the commutative limit, $at * / is reduced to 

detW(f 1 ,f 2 ,...,f N ,f) 



®N*f 



(4.16) 



detW(h,f 2 ,...J N ) ' 

which just coincides with commutative one [23]. In this respect, quasi-determinants are 
fit to this framework of the Wronskian solutions. 
From Eq. (14.131) . we have a more explicit form as 



U2 = d x (jzw'^w;^ 



(4.17) 



where 

W a :=\W(f 1 ,...,f a )\ aa . (4.18) 

The /-reduction condition (L')<_i = or L l = B>i is realized at the level of the soliton 
solutions by taking ot l s = (3 l s or equivalently a s = e(3 s for s — 1, • • • , N, where e is the l-th 
root of unity. 

The present discussion is straightforwardly applicable for NC versions of the matrix 
KP hierarchy [321 M> 125], the toroidal (matrix) GD hierarchy [331 132 EH ES EH [38] 
and the (2-dimensional) Toda lattice hierarchy [39] formulated by pseudo-differential op- 
erators, because on commutative spaces, their exact soliton solutions are described by 
determinants of (generalized) Wronski matrices. 

For example, we can give exact iV-soliton solutions of the NC toroidal ZKdV hierarchy 
(/ > 2)0 which is defined as follows. 

First, we introduce two kind of infinite variables x = (xx,X2, ■ • -) and y :— (yo,yi,y2i, • • •) 
with (x,y) = (xi,yo). Noncommutativity is introduced into these coordinates. Next let 
us define two kind of Lax operators with respect to x, that is, an l-th order differential 
operator P = (L) l >0 and a 0-th order pseudo-differential operator Q, where the coeffi- 
cients depend on the two kind of infinite variables. An differential operator C m i is also 
introduced in terms of P and Q as C m \ := —{P m * Q)>q. Then we can obtain the NC 
toroidal ZKdV hierarchy: 

^-[B m ,Pl, ^-=[B m ,Q-d y l, (4.19) 



dP dO 

= [P m dy + C ml , PI, J^-= [P m dy + C mh Q - dyl . (4.20) 



For I = 2, this includes the NC Calogero-Bogoyavlenskii-Schiff equation [40J. 
The A^-soliton solution is given by 

P = <S> N *d l x $] v \ Q = (d y $ N )*$] f 1 , (4.21) 

where the arguments in $/y is modified as follows: 

f s (x,y) := e^ s ^ + a s e\^ s ^\ (4.22) 
£ r (x,y;a) := £(x; a) + r£(y; a) 

= x\a + x 2 a 2 + x 3 a 3 H h ry + ryi<x + ry 2 ia 21 H , (4.23) 



4 Toroidal IKdV hierarchy is one of generalizations of IKdV hierarchy and first studied by Bogoy- 
avlenskii [33] for I — 2 and developed by Billig, Iohara, Saito, Wakimoto, Ikeda and Takasaki where 
the symmetry of the solution space is revealed to be described in terms of a toroidal Lie algebra, that 
is, a central extension of double loop algebra Q tor = SC\ or [3S1 [3S1 [37] . Hence we call it toroidal IKdV 
hierarchy or toroidal GD hierarchy in the present paper. 
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with ot\ = /3 g , where r is a constant. The proof is the same as the commutative one. (For 
the details, see section 5.1 in |37J.) A key point of the proof is to show the evolution 
equations of <3>at: 



<9$ 



0x m 
d<S> N 



V -($ N d?$J f 1 )<- 1 *$ N = B m *<f> N -Q N dZ, 
(P m * Q)<_! * $ w = (P m d y + C w *) * $at, 



where the following property of quasideterminant plays crucial roles: 

&N*f.= \W(f 1 ,...,f N ,f a )\ N+ i iN+1 = 0, for s = l,---,N. 

This hierarchy generally gives rise to (2 + l)-dimensional integrable equations where space 
and time coordinates are (x, y) and some other coordinate, respectively. 

5 Asymptotic Behavior of the Exact Soliton Solu- 
tions 

In this section, we discuss asymptotic behavior of the multi-soliton solutions at spatial 
infinity or infinitely past and future. In the star-product formalism, all coordinates are 
regarded as c-numbers and hence we can plot the configurations as usual and interpret 
the positions of localized wave packets, effect of coordinate shifts and so on as usual. Here 
we restrict ourselves to NC KdV and KP hierarchies, however, this observation would be 
also applicable to other NC hierarchies. 

First, we present some special properties of the star exponential functions relevant 
to behavior of NC soliton solutions. In this section, we restrict ourselves to a specific 
equation on (2 + 1) or (1 + l)-dimensional space-time and noncommutativity should be 
introduce to some two specific space-time coordinates. Let us suppose that the specified 
NC coordinates are denoted by x% and Xj (i < j) which satisfies = id. 

First, let us comment on an important formula which is relevant to one-soliton solu- 
tions. Defining new coordinates z := x\ + vxj, z := Xi — vxj, we can easily see 

f(z)*g{z) = f(z)g{z) (5.1) 

because the Moyal-product (12.21) is rewritten in terms of (z, z) as [H] 

f(z, z) * g(z, z) = e ™m<d z „-8 z ,8„) f{z ^ z') g ( z », - z ») z , = zll = z (5.2) 
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Hence NC one soliton-solutions are essentially the same as commutative ones for both 
space-time and space-space noncommutativity cases. 

When f(x) is a linear function, the treatment of el^ becomes easy as follows: 

(ef^))- 1 = e ~^ a \ (5.3) 
d x e^ a) = aei^' a) . (5.4) 

The proof can be seen as well from the fact that because of Eq. (15.21) . the star exponential 
function of a linear function itself reduces to commutative one, that is, e* = e^ x,a \ 
These formula are crucial in discussion on asymptotic behavior of iV-soliton solutions. 
Furthermore, the Baker- Campbell-Hausdorff (BCH) formula implies 

e <=(3«) * e «(^) = e (V2)0(a^-c^ e p;a)+£(x;/3) _ = ^(a* -o? &) e tf$\P) ^ e |(^«). (5.5) 

The factor (i/2)9(a l f3 j — (3 l ) can be absorbed by a coordinate shift in a), and hence 
there is a possibility that noncommutativity might affect coordinate shifts by the factor 
such as phase shifts in the asymptotic behavior. When coordinates and fields are treated 
as complex, such a coordinate shift by a complex number causes no problem. However, 
if we want to apply NC integrable equations to real phenomena, such as, shallow water 
waves, then it becomes hard to interpret physically. Let us see what happens in the 
asymptotic region. 

5.1 Asymptotic behavior of NC KdV hierarchy 

First, let us discuss the NC KdV hierarchy and the asymptotic behavior of the iV-soliton 
solutions. The NC KdV hierarchy is the 2-reduction of the NC KP hierarchy and realized 
by putting f3 s = —a s on the iV-soliton solutions of the NC KP hierarchy. Here the con- 
stants a s and a s are non-zero real numbers and a s is positive. Because of the permutation 
property of the columns of quasi-determinants (cf. section 1.1 in |18]). we can assume 
OL\ < a.% < ■ ■ ■ < «AT. 

In the NC KdV hierarchy, the x 2n -th flow becomes trivial and in the x 2n+1 -th flow 
equation, space and time coordinates are specified as (x,t) = (xi,x 2n +i)- 

Now let us define a new coordinate x := x + aj n t comoving with the 7-th soliton and 
take t -> ±00 limit@ Then, because of x + a^t = x + a} n t+{a* n -aj n )t, either e ^ x+a ^ 



5 Such kind of observation for soliton scatterings in NC integrable equations is first seen in [32] ■ (See 
also g].) 
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or e * as( - x+as ^ goes to zero for s ^ I. Hence the behavior of f s becomes at t — > +00: 



fs(x) 



and at t 



-00: 



fs(x) 



- as (x+a 2 a n t) 
a I {x+a 2n t) -ai(x+aj n t) 

e* + die* 

a s (x+a 2n t) 



a a (x+a 2 s "t) 

ar(x+a 2n t) ~aT(x+a 2 "t) 

e* + die* 



a.,e, 



-a„(x+a 2n t) 



S<I 
S = I 
s > I, 

S<I 
S = I 
s> I. 



(5.6) 



(5.7) 



We note that the s-th (s 7^ /) column is proportional to a single exponential function 
e ±a a (x+a s t) ^ Q (|5.4|) . Because of the multiplication property of columns of quasi- 
determinants (cf. section 1.2 in [18]), we can eliminate a common invertible factor from 
the s-th column in \A\ij where s 7^ j. (Note that this exponential function is actually 
invertible as is shown in Eq. (|5.3p .) Hence the iV-soliton solution becomes the following 
simple form where only the J-th column is non-trivial, at t — > +00: 



1 



—OL\ 
(-Ol 

and at t 



\N-1 



1 

-a/_i 



e* + aie* 



a/(eP w) 



N-i n --Zfimh 



1 

oci+i 



„,N-1 

a I+1 



af(ef ( * w) 



+ (-1)^6 



N 



1 

«7V 



„,JV-1 

a N 



a 



./ 

y(jV-l) 



/( A0 



-OO: 



®N*f 
1 

«1 



1 

aj-i 



g(x;o/) — g(x;o/) 



cm e* 



1 

-aj+i 



2V-1 .JV-l 



a 



„N/Mx;°'i) 



+ (-1) 
+ (-1 



AT-1 



iJV-l 



1 

-a at 



/ 



AT I j{N) 



0(i [e* + {— l ) a i e * ) {—oii+i) •-• {—on 
Here we can see that all elements in between the first column and the iV-th column 
commute and depend only on x + aj n t in £(x; ai), which implies that the corresponding 
asymptotic configuration coincides with the commutative oneJ§ that is, the J-th one- 
soliton configuration with some coordinate shift so called the phase shift. The commuta- 
tive discussion has been studied in this way by many authors, and therefore, we conclude 



6 Note that because / is arbitrary, there is no need to consider the products between a column and 
the (N + l)-th column. This observation for asymptotic behavior can be made from Eq. (|4.17p also. 
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that for the NC KdV hierarchy, asymptotic behavior of the multi-soliton solutions is all 
the same as commutative one, and as the results, the N-soliton solutions possess N lo- 
calized energy densities and in the scattering process, they never decay and preserve their 
shapes and velocities of the localized solitary waves. The phase shifts also occur by the 
same degree as commutative ones. 

Finally, we make a brief comment on the 2-soliton solutions. In this situation, space- 
time dependence appears only as two kind of exponential factors e ±a ^( x + a \ nt ) an d e ±Q 2(>+a! n t)_ 
Noncommutativity of them could have effects by the factor e M l / 2 ) a ^ a ^{ a \ n -ai n )e Decause Q f 
the BCH formula. However, if the two kind of parameters satisfy 

-axa^a 2 ^ - a 2 2 n )B = 2irk, (5.8) 

where k is an non-zero integer, then, the effects of noncommutativity perfectly disappear 
at the every stage of calculations and the behavior of the 2-soliton solution perfectly 
coincides with that of commutative one at any time and any location. However the 
condition (15.81) is given specially by hand, and the mathematical and physical meaning of 
this observation is still unknown. 

5.2 Asymptotic behavior of NC KP hierarchy 

Now, let us discuss the NC KP hierarchy and the asymptotic behavior of the iV-soliton 
solutions. The space and time coordinates are (x,y,t) = (xi,X2,x n ) and noncommuta- 
tivity is introduced into some specified two coordinates among x, y and t. The specified 
NC coordinates are also denoted by Xi and Xj with [2^,2^]* = id. Here the constants a s 
and (3 S are non-zero real numbers and the constant a s will be redefined later. 

As we mentioned at the beginning of the present section, one-soliton solutions are 
all the same as commutative ones. However, we have to treat carefully for the NC KP 
hierarchy. From Eq. (I4.17p . naive one-soliton solution can be expressed as follows 

u 2 = d x (d x (ei { *' a) + ael {m ) * (ep Q) + ae^y 1 ) 

= d x ((a l + aP l Aef*-> a 'P ) )*{l + aAef*-> a '' 3) )- 1 ) , (5.9) 

where 

ri(x;a,(3) := x{(3 - a) + y((3 2 - a 2 ) + t((3 n - a n ) 

A := eWC-*?). (5.10) 

We note that the factor A can be absorbed by redefining a coordinate such as x — > 
x + (j3 — a)~ x {il '2)8 (a 1 ft — a 3 f3 l ). The final form of the solution depend only on Xi(f3\ — 
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a]) + Xj(j3j — a J j) for NC coordinates and the Moyal products disappear. Hence there 
becomes no dependence of complex numbers, and the one-soliton solution is the same 
as commutative one in this sense. However now we treat the coordinates as real and it 
would be better to redefine a positive real number a which satisfies a = a A -1 , so that 
fi = e* x '°^ + ael^'^ = (l + ae£ x ' a '^ * el^ x ' a \ in order to avoid such a coordinate shift 
by a complex number. 

This point becomes important for the multi-soliton solutions. The constants a s in 
the iV-soliton solution of the NC KP hierarchy should be replaced with a positive real 
number a s which satisfies a s = fi^A" 1 where A s := e^/ 2 )^"*^ -0 *^, because the iV-soliton 
configuration reduces to a (i-th) one-soliton configuration when we set ot s = (3 S = for 
all s&I). 

Let us define new coordinates comoving with the i-th soliton as follows: 

p := x + airy + afH, q := x + fay + f3?~H. (5.11) 

Then the function £(x,y,t; a s ) can be rewritten in terms of the new coordinates as 
£(p, q, x r \ a s ) = A(a s )p + B(a s )q + C(a s )x r where x r is a specified coordinate among 
x, y and t, and A(a s ), B(a s ) and C(a s ) are real constants depending on aj,/3j and a s . 
For example, in the case of x r = t, we can get from Eq. (15. lip 




1 ( Pip - ajq + aM^~ 2 ~ 
Pi-clA -V + q + K _1 - ft' 1 )* 



(5.12) 



and find 



£ = x + a s y + a™ t 

Pi -a B , a s - ai j n _ x atf^- 2 - a n r 2 ) + a^a?" 1 ~ PT l ) \ . 

Pi -Oj j9/-Q/ I Pi -Oil J 

Here we suppose that C(a s ) ^ C(/3 S ) which corresponds to pure soliton scatterings!^ 

Now let us take x r — > ±oo limit, then, for the same reason as in the NC KdV hierarchy, 
we can see that the asymptotic behavior of / s becomes: 

W) ^ { + aiel^ s = I ( ° 3) 

where A s is some real constant whose value is 1 or a s , and 7 S is a real constant taking a 
value of a s or (3 S . As in the case of the NC KdV hierarchy, the s-th (s ^ /) column is 
proportional to a single exponential function and we can eliminate this factor from the 



7 The condition C(a s ) = C{f} s ) could lead to soliton resonances in commutative case. 
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s-th column. Hence in the asymptotic region x r — > ±00, the iV-soliton solution becomes 
the following simple form where only the 7-th column is non-trivial: 



1 

7i 



1 

J1-1 



7f 



1 

7i 



7/-i 
II -I 



a/e, 

„iV-l„£(£;ar) 



1 

7/+1 



1 

7iv 



/ 



jNA&Pi) 



1 

7/-i 



a j e* 
afei iS ' ai) + ai ^et 

1 + a je ^ w,/3l) 



7/+i 
7/+i 



7 iV-l 

7# 



7f- X 



71 



A? 



7i V -- 1 1 
7f-i 



AT-1 



+ 5j# 



1 

7/+i 



7/+1 



7/+i 



1 

7^ 



7^ 
7# 



./ 

f(N-l) 



fi N) 



Here we can see that all elements between the first column and the iV-th column commute 
and depend only on Xj(/3} — a}) +Xj{f3\ — oPj) for NC coordinates, which implies that the 
corresponding asymptotic configuration coincides with the commutative one. Hence, we 
can also conclude that for the NC KP hierarchy, asymptotic behavior of the multi-soliton 
solutions is all the same as commutative one in the process of pure soliton scatterings, and 
as the results, the N -soliton solutions possess N localized energy densities and in the pure 
scattering process, they never decay and preserve their shapes and velocities of localized 
solitary waves. Asymptotic behavior of two-soliton solution of NC KP equation studied 
by Paniak [T3] actually coincides with our result for n = 3, N = 2. 



Now we restricted ourselves to the NC KP hierarchy, however, this observation would 
be also true of other kind of NC hierarchies, such as, the (2-dimensional) NC Toda lattice 
hierarchy |43j . the NC toroidal GD hierarchy presented in section 4 and the NC matrix 
KP hierarchy [Hj because the soliton solutions could be represented by such kind of 
(generalized) Wronski matrices here, and the asymptotic analysis would be almost the 



same. 



6 Conclusion and Discussion 

In this paper, we studied exact multi-soliton solutions of NC integrable hierarchies, includ- 
ing NC KP and toroidal KP hierarchies and the reductions, in terms of quasi-determinants. 
We found that the asymptotic behavior of them could be all the same as commutative ones 
in the process of (pure) soliton scatterings. This implies that the exact soliton solutions 
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are actually solitons in the sense that the configuration has localized energy densities and 
never decay, and the phase shifts also appear by the same degree as in the commutative 
case. 

It would be reasonable that there is no difference in asymptotic behavior of pure soli- 
ton scatterings on between commutative and NC spaces, because in asymptotic region, 
star-products reduce to ordinary commutative products and the effect of noncommuta- 
tivity disappears. These results imply that we cannot detect effects of noncommutativity 
of space-time by observing such soliton dynamics. However, total behavior of them is 
unknown and it is worth studying further to find different aspects of the NC soliton 
dynamics from commutative ones. 

Dynamics in soliton resonances is also interesting. From the present results of pure 
soliton scatterings, we could naturally expect that the configurations in soliton resonances 
would not be affected by noncommutativity in asymptotic region though we might need 
to make further modifications in the multi-soliton solutions such cLS Qig — (Xg A s 1 for the 
NC KP equation. Quantum treatments of the soliton scatterings is also interesting, such 
as properties of factorized S-matrix of NC sine-Gordon model |45j. Furthermore, the 
existence of multi-soliton solutions is important in integrable systems and the present 
observations might be a hint to reveal NC Hirota's bilinearization, theory of tau-functions 
and the structure of solution spaces. 
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